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Simplified Shear Solution for Determination of the 
Shear Stress Distribution in a Composite Panel from 

the Applied Shear Resultant 

The “simplified shear solution” method is presented for approximating the 
through-thickness shear stress distribution within a composite laminate or panel 
based on laminated beam theory. The method does not consider the solution of a 
particular boundary value problem; rather it requires only knowledge of the 
global shear loading, geometry, and material properties of the laminate or 
panel. It is thus analogous to lamination theory in that ply level stresses can be 
efficiently determined from global load resultants (as determined, for instance, by 
finite element analysis) at a given location in a structure and used to evaluate the 
margin of safety on a ply by ply basis. The simplified shear solution stress 
distribution is zero at free surfaces, continuous at ply boundaries, and integrates 
to the applied shear load. Comparisons to existing theories are made for a variety 
of laminates, and design examples are provided illustrating the use of the method 
for determining through-thickness shear stress margins in several types of 
composite panels and in the context of a finite element structural analysis. 
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HyperSizer automated structural sizing software optimizes 
composite panels based on FEA loads 



• HyperSizer is used extensively by NASA 
for design and trades 

• Employs FEA (NASTRAN/Abaqus) loads 
model to obtain panel and beam level 
loading 

• Uses closed-form panel homogenization 
and localization methods to obtain local 
and ply stresses 

• Evaluates design vs. a multitude of specific 
failure criteria/theories 

- Panel buckling, FS wrinkling, dimpling, core 
crushing, core shear crimping, ~15 strength 

- Considers load cases simultaneously 

• Sizes all geometry of panels/beams and 
selects material for min. weight 

• Updates materials in FEA model to reflect 
new design 

• Iterations between HyperSizer optimization 
and FEA conducted to obtain converged 
design 



Same 


New - HyperFEA automated FEA iterations 
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Previously, through-thickness shear loading was not 
considered for stiffened panels and unstiffened laminates 



Structural FEA 



Candidate Panel Concepts 



I -Stiffened 


£> No TT shear check 


Honeycomb 

Sandwich 


Margin determined 
^ vs. TT shear 
allowable of core 


Solid 

Laminate 


£> No TT shear check 


• In-plane loads (N x , N y , N xy , M x , M y , M xy ) handled consistently 

- Panel homogenization methods and CLT — ► ply-level fields for margin calculation 

• Through-thickness shear loads (Q x , Q y ) only considered for sandwich panel core 

- Cannot determine ply-level margins or perform concept trades 

• Shortcoming identified by NASA - collaborated with Collier R&D to overcome 
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A method was needed to determine ply-level shear 
stresses given the panel-level shear loads 



Desired Characteristics 

• Efficient 

- Will be called thousands of times by HyperSizer during optimization 

• Accurate 

- Integrate to the applied shear resultant 

- Zero at free surfaces 

- Continuous at ply interfaces 

• Able to be applied without knowledge of boundary conditions 

- Not a plate theory solution to a boundary value problem 

- Boundary conditions not available 

• Looking for something analogous to CLT 

- Given laminate layup/ply properties and loads (N, M) 

- Can calculate in-plane ply-level stresses 
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The proposed method is based on the concept of a 
composite beam subjected to bending due to shear 



Z 


vertical cut 

C x 



Actual Transformed 

Cross-Section Cross-Section 



E k 

• Actual ply width (b) replaced by transformed ply width (bn k ), n k =^ 

E 

- Scales the plies by axial moduli ratios 

• Employ plane strain moduli: E k x -+E k x /[i-v k xy v k yx ~] 


• Neutral axis given by: 

TV j TV j 

) A Z^( z t + z k-x ) h n k 

* _ k=l ^ k = 1 ^ 


TV 


TV 

14 


Yjk n k 

k = 1 


k = 1 


Z = 


* 

z — z 


Distance of point from 
neutral axis 
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Beam theory and equilibrium are applied 


• Beam theory axial stress given by, <j x = 

• Consider horizontal cut at z = z. a 


Mz Pxz 


neutral axis 


2n — c 


H 


Equilibrium dictates, r Px % 

-H + ]a x dA = 0 ► H = — \zdA 


Substituting, 


Px "r 

dA-b n k dz ► H = -y- J bn k zdz 


• Integral denoted by Q’, 

z c z c z k Z N -\ 

Q'(k,z) = | bn k zdz = J bn k z dz + J bn k+l zdz + ...+ J bn N zdz 


bn. 


b ni 


Layer k 

A - 



B 



:df 


cross-section 

vertical cut 


Q , {k,z) = dPJz(z]-z 2 k ) + ]- Y, bi 


^ m=k + 1 


- Piece-wise parabolic and continuous at interfaces through the thickness 




horizontal cut 
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Shear stress distribution is determined from shear 
force (H) distribution 



• Shear force distribution, H = —Q'(k,z) 

• Shear stress distribution, *■=— ^-=— e'(*,z) 

xb lb 


• Shear resultant, Q = p/b 


Final form, 



• With, / = X 


k = 1 


12 


bn k t k +bn k t k 


/ A /\ \ 2 


• Applicable to x xz and x yz by switching x and y coordinates 

- No coupling between x xz and x yz 

- Can superimpose effects of Q x and Q y 

• Assumes plate is thin - dominated by bending 
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SSS shear stress fields compared to other 
theories for several layups 



• Compare to Higher-Order Theory for Functionally Graded Materials 
(HOTFGM ) 1 plate solution 

z a x = 0 1 



Symmetry Free T xz = 1 MPa 


i 

1 mm 
I , 

/ / \ 

^ — Pinned / 


— "1 11 ir\r\ ir\r\ - 

\ u mm 


• Compare to HyperSizer Joints 2 bonded laminate in cylindrical bending 
solution 


^boudi, J., Pindera, M.-J., and Arnold, S.M., “Higher-Order Theory for Functionally Graded Materials,” Composites Part B , Vol. 30, 1999, pp. 
777-832. 

2 Zhang, J., Bednarcyk, B.A., Collier, C.S., Yarrington, P.W., Bansal, Y., and Pindera, M.-J., “Analysis Tools for Adhesively Bonded Composite 
Joints Part 2: Unified Analytical Theory,” AIAA Journal , Vol. 44, 2006, pp. 1709-1719. 
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SSS matches well for various 


layups with isotropic 
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SSS matches well for various layups with isotropic 
plies 



Addition of Poisson Mismatch 


0.5 

0.4 

0.3 

0.2 

0.1 

0.0 


- 0.1 


- 0.2 


- 0.3 


- 0.4 


- 0.5 
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Comparison to a true plate theory shows the ability of SSS to 
approximate the shear distribution without knowledge of BCs 



* q<*> 


1 


1 


Compare to Williams’ 3 multilength scale plate theory 

- Solves plate boundary value problem ^ 

- Subjected to sinusoidal pressure loading { 

- Shown to correspond to Pagano’s 4 exact solution 

- SSS applied Q x determined by integrating Williams’ theory x xz distribution from a 
location in the plate (Q x = 0.32 N/mm) 

Material properties: 




06 g - 


X 


Ply Angle (°) 

-30 

0 

30 

45 

60 

90 

E x (GPa) 

2.192 

25 

2.192 

1.325 

1.068 

1 

E y (GPa) 

1.068 

1 

1.068 

1.325 

2.192 

25 

As 

0.4082 

0.25 

0.4082 

0.3245 

0.1989 

0.01 


Layups considered: 

- [0°/90 o /0 o /90 o /0°], [0°/90 o /45 o /90 o /0°] , [-30°/90 o /45 o /60 o /30°] 
-Ply thicknesses: [0.2, 0.4, 0.2, 0.1, 0.1] mm 


^ asymmetric 


3 Williams, T.O., “A Generalized Multilength Scale Nonlinear Composite Plate Theory with Delamination,” International Journal of Solids and 
Structures , Vol. 36, 1999, pp. 3015-3050. 

4 Pagano, N., “Exact Solutions for Composite Laminates in Cylindrical Bending,” Journal of Composite Materials , Vol. 3, 1969, pp. 398-41 1 
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SSS s approximate compared to true plate theory, and 
presence of anisotropic plies causes some deviation 



[07907079070°] 


— Simplified Shear Solution 
-a- HyperSizer Joints 
-•-Williams (1999) 



[0 o /90°/45790 o /0°] 


[-30790745760730°] 


- 0.2 


-0.4 


0.50 



- 0.6 



0.50 


Computational effort of SSS far 
less 

Plate theories not well suited for 
implementation in HyperSizer 
due to presence of BCs 
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Stiffened panel design in HyperSizer is accomplished 
by considering the ply margins due to TT shear 



• Ply-level x xz and x yz are transformed to ply coordinates: 


k 3 l 

cos# 

sin# 


Xtz 

r = 



< 

> 

l r 23j 

-sin# 

cos# 


Im 


TT shear interaction equation considered as failure criterion: 


f 


'13 


Fsu 


+ 


13 J 


'23 


yFSU 23 j 


= 1 




Margin of safety given by: 


MS = 


( 


T 

YJ 

V' 1 oW 13 J 


'13 

Fsu, 


^ 2 f \ 2 
T , 


+ 


'23 

\ FSU 23 V 


1 


- MS = 0 at failure, positive margin good, negative margin indicates failure 


Optimum design is lightest panel that provides all positive margins 

- HyperSizer checks a multitude of margins (e.g., Panel buckling, FS 
wrinkling, dimpling, core crushing, core shear crimping, -15 strength) 
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Determination of TT shear margins in a unstiffened 
laminate 



• Consider [45°/-45 o /0 o /9070790 o /45 o /-4570°] s laminate 

• Ply thickness = 0.005 in. 

• Subjected to Q x = Q y = 200 Ib/in 

• Graphite/Epoxy ply material properties: 


Ej (Msi) 

E 2 (Msi) 

V 12 

Fsu 13 (ksi) 

Fsu 23 (ksi) 

23.35 

1.65 

0.32 

14.8 

5.32 
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Determination of TT shear margins in an unstiffened 
laminate 




Shear Stress (psi) 



Margin positive 
throughout - 
design acceptable 
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Determination of TT shear margins in a honeycomb 
sandwich panel 



Facesheet laminates same as previous unstiffened laminate 
Hexcell HRH 10 (Nomex/Phenolic) honeycomb core thickness = 1.5 in. 
Subjected to Q x = Q y = 1000 Ib/in 
Core material properties: 


E x = E y (psi) 

Fsu xz (ksi) 

Fsu yz (ksi) 

0.0001 

0.225 

0.11 
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Determination of TT shear margins in a honeycomb 
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Determination of TT shear margins in an I -stiffened panel 




• For Q x , stiffener aligned with shear force, (smeared) stiffener webs 
contribute 

• For Q y , stiffener webs do not contribute - shear force applied to 
facesheet only 

• Same facesheet layup, [45/-45/0/90/0/90/45/-45/0] s 

- Ply thicknesses of 0.005 in. 

• Upper flange and web of the I -stiffeners are [0/90/45/-45] s 

- Ply thickness of 0.00375 in. 

• Lower flange layup is [0/90/45/-45/0/0] s 

- Ply thickness of [0.005/0.005/0.005/0.005/0.00375/0.00375] s in. 

• Applied shear loading is Q x = Q y = 200 Ib/in. 
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Determination of TT shear margins in an I -stiffened panel 




— * A 

^ Facesheet 

z 

v 

' Upper Flange 

A 

V 

-« Web 


— *y 

s Lower Flange 


► A 



Shear Stress (psi) 




Margin of Safety 

Margin >1 throughout - design overly conservative 
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Design of unstiffened laminates in a pressure shell 
structure considering TT shear 





14 unstiffened laminate panels within a spacecraft pressure shell: 



[Panel #1019] "SEG Core" 

Panel #1020 "SEG Core" 

Panel #1021 "SEG Core" 

Panel #1022 "SEG Core” 

Panel #1023 "SEG Core" 

Panel #1024 "SEG Core" 

Panel #9000 "Sugar Scoop Solid Laminat" 
Panel #9001 "Sugar Scoop Walls" 

Panel #10013 "SEG Core" 

Panel #10014 "SEG Core" 

Panel #10015 "SEG Core" 

Panel #10016 "SEG Core" 

Panel #10017 "SEG Core" 

Panel #10018 "SEG Core" 




Woven graphite/epoxy ply material properties: 


E! 

(Msi) 

e 2 

(Msi) 

V 12 

G 12 

(Msi) 

FtUj 

(ksi) 

Ftu 2 

(ksi) 

Fsu 17 

(ksi) 

Fsu n 

(ksi) 

Fsu„ 

(ksi) 

10.5 

10.5 

0.036 

0.84 

59 

59 

13 

6.9 

6.9 


• The ply thickness is 0.008 in. 

• HyperSizer is used to determine the number of plies (and thus total 
thickness) needed for each composite panel 

• Candidate laminates consisting of only 0° and 45° plies in a symmetric 
layup were considered 

• 45 load cases (pressure, abort, ascent, reentry, chute loads, landing, etc.) 


www.nasa.gov 22 



National Aeronautics and Space Administration 

Design of unstiffened laminates in a pressure shell 
structure considering TT shear 



Panel thickness 


TT shear failure not checked 



92 lbs 


(in) 


0.16 

0.246 

0.332 

0.418 

0.504 

0.59 

0.676 

0.762 

0.848 

0.934 

1.02 

1.106 

1.192 

1.278 

1.364 

1.45 

1.536 



Margin 



TT shear distribution in panel 10016 


Allowable = 6900 psi 



0.16 _ 

Design with TT shear check 0 . 246 1 

0.332 
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Conclusions 



• A simplified method based on laminated beam theory has been presented for 
determination of the shear stress distribution in laminated plates and panels caused by 
applied shear force resultants 

• The method is straightforward and efficient 

- Does not require the solution of elasticity equations in the presence of the corresponding boundary 
conditions 

- Provides an estimate of the shear stress distribution with no knowledge of the details of the 
structure using only panel cross-section geometry, material properties, and the applied load 

- Analogous to in-plane predictions of CLT 

• The predicted shear stress distribution for laminates are shown to be in excellent 
agreement with HOTFGM and HyperSizer Joints for isotropic plies 

• For anisotropic plies, the agreement between the simplified shear solution, HyperSizer 
Joints, and plate theory is satisfactory given the simplicity and efficiency SSS 

• The method has been applied to predict the biaxial shear distributions in unstiffened 
laminates, a honeycomb sandwich panel, and l-stiffened panel 

- The shear stress distributions in the local ply coordinates are obtained by the relevant 
transformations, enabling determination of the margins of safety of these composite structures 

• Limitations include: 

- Predictions will be less accurate as laminate becomes thick and is less dominated by bending 

- Method predicts zero shear stress throughout laminate if applied shear load is zero 

- No interaction between x xz and x yz - determined independently and superimposed 
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